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Abstract 



The effective equation of motion that describes the different monomer alter- 
nation along the heteropolymer chain is proposed. On this basis the super- 
symmetry field scheme is built up to analyze memory and ergodicity breaking 
effects. 
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Heteropolymer represents an ensemble of chains of different type monomers (many-letter 
sequences), that can experience the freezing and the microphase separation (MS) (see [1] 
and references therein). Firstly, MS was studied within random phase approximation [2] 
for the case of block copolymer AB with arbitrary fraction / of monomers of type A. At 
/ 7^ 0.5 with temperature decrease the system undergoes MS sequence of the first order, 
where period 2n/ko of space structure is of order of block length and does not depend on 
temperature. By virtue of condition ko ^ the order parameter fluctuations give divergent 
contribution to thermodynamic values, so that continuous phase transition transforms to 
the weak first order transition [3]. Passing to random heteropolymers both MS and freezing 
transitions remain, but here the space period is strongly dependent on temperature [1, 4]. 
According to the field considerations [5, 6] the fluctuations suppress both transitions in 
random copolymers, but copolymer melts in disordered media suffer these transformations. 

Similar to spin glass [7] , the main tool of statistical mechanics of random heteropolymers 
is the replica trick. Besides, the transfer matrix method [8] and kinetic approach [9] have 
been used. At the same time, it is well known in theory of spin glass that within the frame- 
work of Sherrington-Kirkpatrick model the replica method is identical to the supersymmetry 
(SUSY) field approach [10,11]. Being based on introducing Grassmann variables, the SUSY 
field components have explicit physical meaning, so that in contrast to the replica trick the 
SUSY method arrives at physically interpreted results [12]. Up to now, the SUSY method 
in polymer theory has been proposed by Vilgis [13] and did not obtain further development. 

Usually, the SUSY method is based on the equation of motion of Langevin type. But 
covalent bonds, that ensure the formation of polymer chain, make such way ineffective be- 
cause the dynamic theory of polymers is much more complicated than statistical mechanics 
of usual many-body systems [14]. Therefore, one needs the effective motion equation in- 
stead of dynamical one. The main aim of this work is to obtain such equation for random 
heteropolymers and to build up corresponding statistical SUSY scheme. 

Let us start with the directed polymer that represents Gaussian chain with the proba- 
bility \I/(R, A) to find the end-point A at the point R. As is known, the function ^(R, A) 
obeys the Schroedinger-type equation with imaginary time — iA [14] 

dm/dN = (Dd 2 /dB? - U{R, A)) (1) 

where number A ^> 1, D is the effective diffusion coefficient, C/(R, A) is the external field. 
In the limit where A — > oo the solution of Eq.(l) can be written in the form of functional 
integral over dependence r(n) of chain coordinate on number of internal monomers: 

tf(R, A) = J exp(- 1 S RiV {r(n)}/2 J D)5r(n). (2) 

Here the action S(TL, A) = SRjv{r(n)} = L (r(n))dn with the fixed end-points r(0) = 
and r(A) = R is determined by Lagrangian of the Euclidean field theory [10] 

^(^r)W»), < 3 > 

where the effective kinetic energy, for which the continuum limit r(n+l)— r(n) — > dr(n)/dn is 
assumed, is caused by the covalent bond between monomers of polymer chain [14] . Inserting 
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Eq.(|) to Eq.(l) gives the Jacobi-type equation OS/dN = Dd 2 S/dR 2 - (1/2) (OS/dR) 2 + 
2DU . After introducing the generalized momentum p = dS/dR and the total deriva- 
tive dp/diV = dp/dN + (pd/dR)p the last equation takes the linear form by Burgers: 
dp/diV = D (d 2 p/dR 2 + 2dU/dR). The above relations are well-known in theories of di- 
rected polymers, stochastic growth and kinetic roughening phenomena (see [15]). 

Our basic observation [16] is that the Schroedinger-type equation ([!]) takes the form of 
the Fokker-Planck equation 

2 ^ J 



dN \ OR 2 

for the probability distribution 

P(R, N) = (R, N) exp{-V(R)/2D}, (5) 

where the dependence on R is determined by the effective potential V = — / FdR. 
Corresponding force F is related to the initial potential U in Eq.(|l|) as follows: U = 
F 2 /4D + (l/2)d¥/dR. As is well known [17], Eq.(f|) determines the probability to real- 
ize solution of corresponding Langevin-type equation 

dR/dN = F(R, N) + ((N) (6) 

for stochastic function R = R(N). Here the stochastic force ( is subjected to the white-noise 
conditions (((N)) = 0, (((N)((N')) = 2D5(N - N'), where the angular brackets denote 
averaging with respect to the distribution (0). 

In going from the case of directed polymer to the random heteropolymer AB, the coordi- 
nate R of the end-point iV becomes a stochastic Ising variable 9(n), where 8(n) = 1 if n-th 
segment is of type A and 9(n) = — 1 otherwise. Corresponding to the quenched disorder 
in fixed sequence of different type segments, the law 9{n) of monomer alternation along 
the chain is described by the master equation analogous to that determining the Glauber 
dynamics [|T^]. Relevant sequence correlator a(n)a(n') for effective spin u{n) = 6(ri) —9(n), 
being deviation of microscopic value 9(n) from the average 9(n), takes the form || 



a(n)a(n') = C 2 exp(-\n - n'\/l), (7) 
C 2 = 4/(1 - /), f = (l/2)(l + 9(nj), 

where overbar denotes the averaging over composition sequence (quenched disorder), I is the 
correlation length and / is the fraction of type-A monomers. 

Possessing the exponential correlator (|7|), stochastic variable a{n) is governed by the 
effective motion equation da/dn = —a/l + s(n) with the white noise: s(n) = 0, s(n)s(n') = 
2C2l~ 1 5(n — n'). Relation between the microscopic value cr(n) and stochastic 5-correlated 
variable s(n) is given by the equality [17] o{n) = f$ e~( n ~ m ^ l s(m)dm. Contrary to the 



colored noise cr(n), the white noise s(n) has the Gaussian distribution function determining 
quenched disorder with the intensity Cy -1 . Then the locally averaging field 

T](r,n) = {4C 2 )- 1/2 a(n)5(r - r(n)) (8) 

represents the order parameter (hereafter the monomer volume is assumed to be equal unity). 
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In addition to the terms of type included in in Eq.@, the effective equation of motion for 
field (|) must contain the inhomogeneity contribution Dd 2 r]/dr 2 that has Fourier transform 
—Dk 2 ^. As a result, the equation of motion for the Fourier transformation ??k(n) of field 
(§) takes the Langevin form: 

dr) k /dn = -Dk 2 ^ - &H/dr]* k + ( k . (9) 

Here, as above, the continuum approximation is used for effective time n 3> 1, and the 
definition of effective force / k = —dH/dr/^ is taken into account. The white noise Ck = Ck(^) 
is defined by the conditions 

(Ck) = 0, (Ci(n)Cf(n')) = 2<W(^ - ri), (10) 

where the angular brackets denote averaging over chain conformations (thermal disorder). 
Taking into account the effect of fluctuations, the effective Hamiltonian reads |l],|5]|| 

W = X^kl^kl 2 - 7^ w kk'|?7k| 2 |77k'| 2 + fv (r)dr; (11) 
k 1 kk' J 

r k = r + 2D(k - k ) 2 } k^ = 2{itD) l ' 2 l(2r) l l\ 
r = r+ (3/4vr)/- 2 (2r)- 1 / 2 , r = T 1 - C 2X ; 
Wkk , = 4a 2 /^ 2 (iV J D)- 1 (k 2 + k' 2 )- 1 ; 
v = -(/i/3!)^ 3 + (A/4!),? 4 , 

fi = 12C 3 C 2 1/2 r\ A = 24(1 + 5C%/C 2 )r 1 , 
C 2 ee 4/(1-/), C 3 =|l-2/|. 

Here a, x are the interreplica overlapping and the Flory parameters, correspondingly. 

The following application of the SUSY scheme is straightforward ||10|| . One has to intro- 
duce the generating functional 



SH/Srfc ee dH/dr)* k + 2D(k - fc ) V 



J det 


^Ck 







(12) 



being the average over noise (k(n), where 5-function accounts for the motion equation 
(9), the determinant is Jacobian of transformation from £ k to Then the functional 
Laplace representation is used for ^-function, that gives a ghost field y? k (n). One needs 
to use Grassmann conjugated fields ipk(n), ^(n) |L0| in order to write the determinant in 
exponential form. Then, assuming that conformation averaging in Eq.(|T2"D is Gaussian with 
variance 1 (see Eqs.flTO])), the standard form is derived (cf. Eqs.(2,3)) 



Z{ v } = J P{t], <p; tfj, tp}5^5 2 tfj; 5 2 tp = 5^5tp, (13) 

/•N 

P{V^;^^} = exp(-^{r/,y?;^,^}), S = Ldn, 

Jo 

[(<pf) + &j) - V 2 /2) + (H'{r)}(p + i)H" {r,}^)] dr. 
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Here, the point denotes the derivative with respect to "time" n, the prime denotes the 
functional derivative with respect to the field (8). 

The last expression in Eqs.flllf) takes the simplest form 

L = ~ J A($)d¥, A = $>;;PZ>$ k + ft{$ k }, (14) 

k 

V = d/d$ + <&8/dn, V = 8/dd + dd/dn 

within the SUSY field representation 

$ = r] + ip<d + dip + (15) 

where Grassmann coordinates i?, d obey the same relations as fields ip, ip. Here the new 
field 0k = 77k — V^k is introduced, the functional 7i{^} has the same form as the effective 
Hamiltonian (|TTD, where the order parameter 77k is replaced by the SUSY field $k, Eq.(^5|), 
and the term (5TC/5r)l)rjk + c.c. is the total derivative d7i{7?k(7i)}/dn that can be omitted. 
According to []T2], the physical meaning of the SUSY field components is as follows: <p> = 
—STC/5r]* is the field being conjugated to the order parameter 77, ip = 77 — (p is the most 
probable value of fluctuation of the conjugate field, and the combination ipip gives the density 
of sharp interphases. So, the using of the whole 4-component SUSY field flT5[) corresponds to 
the strong segregation limit. In what follows we concentrate on the simple case of the weak 



segregation limit where ipip = 0. Then, the SUSY field (15) is reduced to the 2-component 
form 

$ = 7i + 6(p, (16) 

where self-conjugated nilpotent variable 9 = dd is introduced. As a result, the Lagrangian 
(ITJ) takes the form 



L = i / A($)d#, A = ^$ k D$ k + H{$k}; (17) 

D = -d/dd + (1 - 268/ 86) d/dn. 

Here the fluctuation field (p is used instead of conjugate field and the lowest power of 
derivative 8/8n is retained. To correspond to Lagrangian (|T7|) the motion equation for 
SUSY field (USD reads 



D$ k = -5H/5&1. (18) 

Let us introduce now the SUSY correlator 

C k (n, 6- n', 6') = (%{n, 6)$ k (n', 6')) #(n - n'), (19) 

where $(n — n') = 1 if n > n' and d(n — n') = otherwise. Multiplying Eq . (|i~8|) by $ k and 
averaging the result, one gets within zeros approximation (w = v = in Eqs.(|TT|)) (TJ 



(20) 

r k + v 
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Here, conventional frequency v denotes Fourier transformation over "time" being the 
monomer number n. The most important feature of expression fl20|) is the characteric struc- 
ture with respect to combination of the nilpotent variables 9, 9', that is inherent not only 
zeros approximation but arbitrary SUSY correlator. To this end it is convenient to introduce 
basis supervectors 

A(9, 9') = 9, B{9, 9') = 9', T(9, 9') = 1, (21) 

containing functional product X(9, 9') = J Y(9, 9") Z(9", 9)d9" for any SUSY vectors X, Y, 
Z. It is easy to see, that basis SUSY vectors (|2l"D obey the following multiplication rules: 
A 2 = A, B 2 = B, BT = T, TA = T, other products are zero. Since the set of vectors A, 
B, T is closed, it is convenient to expand any SUSY correlator over this basis: 

C = G+A + GB + ST. (22) 

Hereafter the subscripts k, v are suppressed for brevity. Using Eqs.(|i~6"D, ([19]), one gets for 
coefficients of expansion Q2"2]): G + = ((f*r)), G_ = (rfip), S = (\r]\ 2 ). So, G± represent 
advanced and retarded Green functions and S is the structure factor. In accordance with 
Egs-POj), pi|), (^2]), within zeroes approximation these functions are 

Gf = (r ± iu)-\ = GfG^ = (r 2 + u 2 )' 1 . (23) 



The Dyson equation for the SUSY correlator (19) reads [12 



CT 1 = (C^-wC-S, (24) 

where w = 2 ■ Q 1 ^n^ 4 ^D^ 1 {a /l) 2 is the typical value of the sum YJ k , w^/ in Hamiltonian 
(PI), S is the SUSY self-energy operator. Using Eqs.(p3|) and expansions (p^), 

S = S + A + S_B + ST, (25) 

the Dyson equations for the components G±, S become 

G ± x + wG± = (r ±h/) - E ± , (26a) 

S = (1 + 2ttC 2 1- 1 5{u) + S)G + G_(1 - wG+G-Y 1 . (26b) 

Here 5-term appears to take into account the contribution to Lagrangian due to quenched 
disorder. Coefficients of SUSY expansion ( pB"D are [pl|] 

s±(n) = (V + Y s ( n n 5 H G ±H' ( 27a ) 

E(n)= L 2 + ^S(n)\s 2 (n), (27b) 
where the r-representation is used for macroscopically homogeneous system. 



6 



Let us introduce now the memory parameter q = (r](n = N)r](n = 0)) and the non- 
ergodicity parameter A = g — g that is the difference between isothermal susceptibility go = 
G-{y = 0) and thermodynamic value g = G_(z/ — > 0). Being coefficients in SUSY expansion 
(^2|), the main correlators of non-ergodic system with memory acquire the elongated form: 

G ± (i/) = A + G ±0 (i/), S(n)=q + S (n), (28) 



where the index denotes the components corresponding to ergodic system without memory. 
Inserting Eqs.(^) to Eqs.(27) provides the self-energy components T,±(n), £(n) as functions 
of parameters A, q within the "time" -representation. Since the Dyson equations (26) require 
"frequency" Fourier transforms, it is convenient to use the fluctuation-dissipation theorem 
So(n — > 0) = G±o( z/ - ► 0) = g, T,±o(v — >• 0) = T,o(n — » 0). As a result, the Dyson equation 
(26) gives in the //-representation 



1 



2 1/2 \ 2 

w % ~2y + Y q j qg ° 



C2l~ l gl 



(29) 



c (1 + S )G + G- 

^ l-[ w + + A 2 g/2)g]G + G_ • 1 Ja) 

The first of these equations corresponds to the 5-term due to memory effects, the second 
one corresponds to "frequency" v ^ 0. At v — ► the characteristic product is G + G_ — ► g 2 , 
so that the pole of the structure factor (29a) 

w + (V + y<Zj Q = 9o 2 (30) 

determines the point of ergodicity breaking. Analogously, Eq. ( |26a| ) gives the equation for 
the susceptibility g = G-{y — > 0) 



wg 2 + ^rg {g + qf-q 1 +'—g {g + qf-q* =rg-l. (31) 



2 . . A 2 



2 
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The system of equations (29)— completely describes thermodynamic behaviour of ran- 
dom heteropolymer by analogy with the spin glass theory @. So, Eqs.(29), ([H]) play a role 
of Sherrington-Kirkpatrick equations, and Eq. ( |30"D defines the point of de Almeida-Thouless 
instability. 

According to Eqs.(29), fl3D|), the memory parameter is given by the cubic equation 

(lx 2 /2 + \ 2 q/3)q 2 = C 2 r 1 . (32) 

Taking into account definitions ([0]), it is seen that for copolymers close to symmetric com- 
position / = 0.5 (C3 <C C-i) the first term in brackets of Eq . (|3~2"D is negligible, and the 
dependence q oc Z 1 / 3 takes place. In the opposite case of dilute copolymer, where / < 1 
(G 2 <C G 3 ), we have the less value q oc fl 1 ^ 2 . Naturally, the memory parameter diminishes 
always with correlation length decreasing. 

The behaviour of isothermal, g , and thermodynamic, g, susceptibilities is described by 
the equations (|30D and (0), respectively. The magnitude g depends only on the quenched 
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disorder I, whereas the latter, g, is defined by both values I, x of quenched and thermal disor- 
ders (see Eqs. flTTD for dependence r(x))- The point of ergodicity breaking, r , is determined 
by equation 

2 \ 2 

2 - rg + ^-gl + —g 3 (g + 3g) = 0, (33) 

that follows from Eqs.(|30D, (PH) at g — go. At the point of MS, r c , one has dg/dr = — oo, 
and Eq.([H]) gives the condition 

A 2 

w + fx 2 (g + q) + —{g + qf = g~ 2 . (34) 
The form of phase diagram given by Eqs.(|33|), (|3^) is shown in Fig.l. 
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FIGURES 

FIG. 1. The phase diagram for random copolymer with a = (a) and <j = 5(b)(Z = l,D = l; 
solid and broken curves are the MS and ergodicity breaking lines correspondingly) 
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This figure "fig.gif" is available in "gif" format from: 



http://arXiv.org/ps/cond-mat/9807063vl 



